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Abstract 

Naturalness in supersymmetry (SUSY) is under siege by increasingly stringent LHC 
constraints, but natural electroweak symmetry breaking still remains the most powerful 
motivation for superpartner masses within experimental reach. If naturalness is the wrong 
criterion then what determines the mass scale of the superpartners? We motivate su- 
persymmetry by (1) gauge coupling unification, (2) dark matter, and (3) precision b — r 
Yukawa unification. We show that for an LSP that is a bino-Higgsino admixture, these 
three requirements lead to an upper-bound on the stop and sbottom masses in the several 
TeV regime because the threshold correction to the bottom mass at the superpartner scale 
is required to have a particular size. For tan f3 « 50, which is needed for t — b — r unification, 
the stops must be lighter than 2.8 TeV when At has the opposite sign of the gluino mass, as 
is favored by renormalization group scaling. For lower values of tan/3, the top and bottom 
squarks must be even lighter. Yukawa unification plus dark matter implies that superpart- 
ners are likely in reach of the LHC, after the upgrade to 14 (or 13) TeV, independent of 
any considerations of naturalness. We present a model-independent, bottom- up analysis of 
the SUSY parameter space that is simultaneously consistent with Yukawa unification and 
the hint for = 125 GeV. We study the flavor and dark matter phenomenology that 
accompanies this Yukawa unification. A large portion of the parameter space predicts that 
the branching fraction for B s — > /i + fj,~ will be observed to be significantly lower than the 
SM value. 
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1 Introduction and Conclusions 



For thirty years Weak Scale Supersymmetry has been a leading candidate for a deeper theory 
underlying the Standard Model (SM); it is motivated by a natural electroweak scale, precision 
unification of the gauge couplings, and Lightest Super Partner (LSP) dark matter from thermal 
freeze-out. However, experiment has called into question the naturalness of the minimal theory, 
first by the LEP2 bound on the Higgs mass and second by direct superpartner searches and 
preliminary evidence for a Higgs boson with a mass near 125 GeV [H [2] from 2011 LHC data. 
If confirmed, the latter implies fine-tuning of the weak scale in the Minimal Supersymmetric 
Standard Model (MSSM) of at least 1% [3J. Furthermore, LSP dark matter from thermal freeze- 
out requires a heavy gravitino, further increasing the fine-tuning. 

Thus LHC data suggests two very different schemes for supersymmetry with precise gauge 
coupling unification: 

• A natural weak scale. An extension of the minimal model is required, for example by 
adding a gauge singlet superfield and/or violating R-parity. LSP dark matter from thermal 
freeze-out may be lost or may require a considerable modification of the theory. 

• A fine-tuned weak scale. The minimal model survives, allowing LSP dark matter from 
thermal freeze-out. 

While both schemes were studied before the LHC, and while there are many avenues for su- 
persymmetry with a natural weak scale, the possibility of a fine-tuned weak scale is growing in 
importance. 

With a fine-tuned weak scale, the supersymmetry breaking scale in becomes decoupled from 
the weak scale, so that discovery of superpartners at the LHC does not appear to be necessary. 
For example, in Split Supersymmetry j3] the scalar superpartners can have a mass m many 
orders of magnitude larger than the weak scale, while only the fermionic superpartners are in 
the TeV domain. Similarly in Anomaly Mediation [5], E] without sequestering, and in Spread 
Supersymmetry [7] , the scalar superpartner mass scale m is larger than the weak scale by powers 
of an inverse loop factor (for early work along these lines see [HI IH])- In all these theories 
supersymmetry may well be out of reach at the LHC. 

Are there any experimental constraints that limit m? A Higgs mass near 125 GeV provides 
a constraint on m, but unfortunately only a weak one. For example, in Split Supersymmetry 
m < 10 5 TeV HDHTJ]. More generally, at one extreme, with the Higgs quartic boosted by large 
top squark mixing, m ~ TeV is possible. At the other extreme, rh ~ 10 13 TeV is possible if the 
Higgs quartic arises entirely from SM loop effects |12| . although this may require a threshold 
to the quartic coupling at the unified scale. Flavor-changing and CP violating effects provide 
similarly weak constraints: if fa > 10 3 TeV, virtual superpartners do not upset the successful 
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CKM understanding of the SM even if the squark and slepton mass matrices have arbitrary 
off-diagonal entries. However, this is clearly not a firm limit since the flavor matrices may have 
a hierarchical structure with small off-diagonal entries, as in the quark mass matrices. 

In this paper we constrain m by requiring precise unification of the b and r Yukawa cou- 
plings at the mass scale of gauge coupling unification [131 HH EE]- Thus, our motivations for 
supersymmetry are the conventional ones, except naturalness is replaced by Yukawa unification: 

• Precise unification of the gauge couplings. 

• Precise unification of the b and r Yukawa couplings. 

• Lightest Super Partner (LSP) dark matter from thermal freeze-out. 

It seems significant that in the MSSM both gauge and Yukawa unification can result from 
very simple unified theories, for example, an 5*0(10) theory with each generation in a 16-plet 
spinor, ip, and the third generation masses arising from the single interaction ip^if)^. If unifies 
the two Higgs doublets of the MSSM, H u ^, then t — b — r unification [161 HZ] occurs, while if <p 
contains unequal components of H u and Hd then b — r unification results. The smaller Yukawa 
couplings involving lighter generations may have a more complicated origin, for example from 
higher dimensional interactions, radiative effects, or environmental selection. 

Ignoring supersymmetric threshold corrections, two-loop renormalization group scaling in 
the MSSM yields b — t unification with a precision of about 20% at tan (3 < 10, and t — b — r 
unification with a precision of about 10% at tan ~ 50. Given that the quark and charged lepton 
masses vary over six orders of magnitude, this Yukawa unification is very striking, perhaps the 
best hint we have of any symmetry structure underlying quark and lepton masses, although 
evaluating its significance is not easy. In figure [TJ the d/e, s/fi and b/r Yukawa ratios are shown 
as a function of scale in the MSSM, with superpartners at 2 TeV and tan (3 = 50. The dashed 
lines ignore the finite (i.e. non-log) supersymmetric threshold corrections, and the shaded bands 
include ±2cx uncertainties on the quark masses. The d/e and s/fi Yukawa ratios at the unified 
scale are close to the Georgi-Jarlskog |20j values of 3 and 1/3, but in this paper we focus on 
the third generation. The solid line shows the effect of including a 12% finite supersymmetric 
threshold correction for the b Yukawa coupling, as needed for precise Yukawa unification. 

We will conclude that the above three motivations for supersymmetry lead tom~ 1 — 10 TeV. 
Yet, once naturalness is ignored (or multiverse arguments are made for the weak scale and the 
dark matter abundance) one might question whether such a supersymmetric theory is preferred 
over the SM with axionic dark matter. It is well-known that gauge coupling unification occurs 
more accurately in the supersymmetric case; here we stress that b/r unification in the SM 
requires a unified threshold correction of 60% - a factor 5 larger than the supersymmetric 
threshold correction required for t — b — r unification. The combination of gauge and Yukawa 
unification is a powerful argument for (multi-) TeV scale supersymmetry. 



3 



Yukawa ratios 



10.0 




m susy = 2 TeV 



100 



10- 



10 8 10 11 
fi [GeV] 



10 



14 



10 



17 



Figure 1: Yukawa couplings ratios for d/e, s/fi and b/r as a function of scale, \i. The three 
dashed lines result from using the central values of the masses derived from experiment, two-loop 
running in the MSSM [18J with all superpartners at 2 TeV and tan /3 = 50, but ignoring the 
finite superpartner threshold corrections. The colored bands allow for ±2a uncertainties in the 
quark masses. The solid curve includes a 12% finite correction at the TeV scale to the b Yukawa 
coupling. The running MS masses evaluated at rriz are taken from Ref. [T9] . 

The key question becomes the origin of the 10 - 20% correction necessary for Yukawa uni- 
fication. Since unified threshold corrections are small for gauge couplings, < 1%, they are also 
expected to be small for Yukawa couplings. Although there is large SU(5) breaking in the s/fi 
Yukawa ratio, this typically affects the third generation b/r boundary condition at the level of 
less than 1%. Thus, while it is possible to construct theories which have the required b/r correc- 
tion arising entirely from either boundary condition or threshold corrections at the unified scale, 
in a very wide class of theories the correction must arise from the supersymmetric threshold. 

A large 10 - 20% correction from the supersymmetric threshold could come from logarithmic 
terms, with at least some superpartners far above the weak scale, or it could come from finite 
threshold corrections. With m at the TeV scale, figure [T] shows that a finite threshold correction 
is needed. We find that raising rh above the TeV scale does not eliminate the requirement of 
a finite threshold correction, whether all superpartners are raised together or Higgsinos and/or 
gauginos are kept at the TeV scale. Furthermore, such large finite corrections are easy to obtain 
since they are enhanced by tan/3, the ratio of electroweak VEVs [IT] . In particular, the leading 



4 



M 3 



g 



a; 



« X > 



H. 



d 



Figure 2: The leading, tan f3 enhanced contributions to the finite part of the bottom Yukawa 
threshold. 

contribution to the finite bottom threshold in the heavy squark limit is 



St = -^^tan/3 - -^^tan/3. (1.1) 
b 12tt 2 ml 32tt 2 m? K ' 

b t 

Two schemes allowing Yukawa unification at large tan /3 were introduced in Ref. [32] . For light 
squarks of a few hundred GeV individual contributions to S( m are typically too large, but there 
could be a cancellation. With somewhat heavier squarks the contributions become suppressed, 
with 5[ m oc fi/m where /i is the Higgsino mass parameter, and can give the desired 10-20% 
correction without a cancellation. Motivated by naturalness, many studies of Yukawa unification 
have focussed on the light squark case (for example Refs. [HI [22j [23J, 121]); although Tobe and 
Wells were early advocates of easing the naturalness requirement, studying Yukawa unification 
with heavier squarks [25] . 

We concentrate on a Higgsino/bino LSP in the MSSM, so that the dark matter abundance 
constrains the Higgsino mass parameter, \i < 1 TeV. As the scale of supersymmetry breaking is 
increased above /x, the finite supersymmetric threshold corrections to b/r become suppressed, as 
they are proportional to \ijm. Since a large 10 - 20% finite supersymmetric threshold correction 
is necessary for exact Yukawa unification, there is an upper bound on m. For tan = 50 we 
will see that this limits the third generation squark mass to be less than about 2.5 TeV when A t 
has the opposite sign of Ms, as is favored by Renormalization Group (RG) scaling As tan/3 
is reduced this limit becomes stronger. Thus the simplest supersymmetric scheme for Yukawa 
unification and LSP dark matter imposes a powerful constraint on m, even when naturalness is 
ignored. When A t and M 3 have the same sign, the limit weakens to about 9 TeV for tan (3 = 50. 



1 The gluino mass, M 3 , appears in the RG for A t (see equation 4.3) and drives it to have the opposite sign 
unless A t is chosen to have a large value in the UV. 



5 



Furthermore, with m < 1 TeV a variety of constraints make the theory implausible. Yukawa 
unification with large tan/3 leads to excessive contributions to b — > 57 and to <5^ m , requiring 
cancellations in both cases. A Higgs mass near 125 GeV is hard to achieve, and LHC bounds 
on squark production become stringent. Hence we are led to Yukawa unification with squark 
masses in the several TeV range. 

Since we are dropping naturalness, one might wonder if there is a supersymmetric spectrum 
having squark masses many orders of magnitude above the weak scale that has both gauge and 
Yukawa unification with small unified thresholds and LSP dark matter. In split supersymmetry 
this is not the CclS6 ; clS the additional logarithmic supersymmetric threshold corrections to b — r 
unification have the wrong sign, while the finite corrections essentially vanish [26]. Having only 
the Higgsinos at the weak scale leads to the same situation, while keeping only the gauginos at 
the weak scale significantly degrades gauge coupling unification. 

The organization of this paper is as follows. In section [2j we perform a model independent 
analysis of b— r unification. As a function of the supersymmetric threshold corrections to the t, b, 
and gauge couplings, we find the parameter choices that give precision Yukawa unification with 
a small GUT-scale threshold correction. In section [3] we discuss various aspects of the GUT-scale 
boundary condition and threshold corrections. We begin in subsection 3.1 by showing that it is 
possible for the first two generations to have non-unified boundary conditions, and to have the 
necessary CKM mixing with the third generation, without generating large corrections to b/r. 
Typically GUT threshold corrections are small, but in subsection 3^ we discuss two possible 
sources of large corrections: (1) a large value of yt near the GUT-scale (which results from low 
tan /3) and, (2) the threshold corrections that arise from Yukawa unification in an orbifold GUT. 

Our detailed exploration of the MSSM parameter space begins in section |4| where we deter- 
mine which parameters yield both precision b — t unification and a 125 GeV Higgs mass. By 
imposing that bino-Higgsino dark matter does not overdose, \i < 1 TeV, we arrive at upper 
bounds on the stop and sbottom masses in the several TeV regime. Having identified the al- 
lowed parameter space, we move on to study its phenomenology in section |5j In subsection 5T 



we consider SUSY contributions to the rare 5-meson decays, b — > and B s — > n + ji~ . For 
the parameters that give Yukawa unification, b — > s'j is an important constraint, but a large 
portion of the parameter space is allowed without requiring fine-tuning. Meanwhile, we find 
that B s — > is observably depleted in much of the parameter space, and therefore provides 



a promising discovery handle. We present the dark matter phenomenology in subsection 5.2 



We determine the relic abundance and direct detection cross-section within the parameter space 
that gives b — r unification. Although the present limits from direct detection are weak, we find 
that XENON1T will probe a sizable fraction of the parameter space. In appendix [A] we consider 
how extra charged states modify Yukawa unification. In appendix [Bj we describe our notation 
for various flavor bases that are useful for understanding rare 5-meson decays. In appendix [O 
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we present an updated, model- independent, study of the parameter space of the well-tempered 
neutralino [27], including recent lattice values for the strange quark content of the nucleon [28J. 



2 Model-Independent Analysis of Threshold Corrections 

In any scheme that embeds the SM in a supersymmetric theory at scale rh and has quark-lepton 
and gauge coupling unification at scale Mu, the couplings will receive threshold corrections at 
both the supersymmetric and unified scales. In this section we present a general analysis of the 
necessary size of the threshold corrections, without appealing to specific supersymmetric spectra. 
Armed with the model-independent results of this section, we will consider explicit spectra in 
section HI 

We match the gauge couplings, g a , with a = 3,2,1 for 577(3), 577(2), U(l) and the third 
generation Yukawa couplings, yi with i = t, b, r of the SM to those of the MSSM by introducing 
threshold corrections 5^ a 

y MSSM {Mz) = vTiMzl (1 + y^ SM (M z ) = y ll^l (1 + 6br) (21) 

sin l) cos jj 

9™ SSM (M z )=g s a M (M z )(l + 5 a ). (2.2) 

Although we are interested in rh up to three orders of magnitude above the weak scale, for 
simplicity we define the above thresholds at the reference scale M z . The unified couplings result 
from matching to the MSSM couplings at Mu 

y(Mu) = yt ISSM (Mu)(l + e)= y™ SSM (Mu) (2.3) 

and 



g 2 (Mu) gl(Mu) 

where 1/g 2 is the average of l/g%. The unified scale, Mu, is defined by minimizing e 2 g = J2 a e a- 
Between the matching scales we evolve the couplings using two-loop MSSM renormalization 
group equations [18] , so that b — r unification leads to a constraint between the thresholds 

e = e(5i, a , tan /3). (2.5) 

When specializing to t — b — r unification, tan (5 is fixed by setting y t = y b at Mu- 

The SUSY threshold corrections, 8 ita , depend on the superpartner spectrum and at one-loop 
are the sum of terms proportional to the log of superpartner masses, 8°%, and terms that are 
independent of logarithms of superpartner masses, 6{™ '. Setting all superpartner masses equal, 
gauge coupling unification is very precise, with e g = 0.013(0.017) for rh = 0.1(100) TeV. For 



7 



superpartner masses in the range of interest to us, gauge coupling unification is remarkably 
insensitive to the supersymmetric thresholds. On the other hand, for the example shown in 
figure [IJ with all superpartners at m = 2 TeV, tan/3 = 50, and 5[ % h n Ta all set to zero, b — r 
unification requires large unified thresholds, e = 0.23. We now analyze more generally the 
constraint on the thresholds, (2.5), required for b — r and t — b — r unification. 




Figure 3: Contours of the GUT-scale threshold e, as a function of the yb and g% supersymmetric 
thresholds, for tan/3 = 10,40. We define e by + e) = y T at the GUT scale, and the red 
region with |e| < 0.02 has precision Yukawa unification. The blue contour shows the logarithmic 
contribution to 5b for a particular spectrum with degenerate scalar masses, described in the 
text. The purple arrow represents the 2a experimental uncertainty on the bottom mass. The 
top Yukawa threshold, 5 t , is chosen to vary with 5^ , also assuming a spectrum with degenerate 
scalars. For simplicity, we have set 5\ = 52 = since these threshold corrections are numerically 
unimportant for Yukawa unification, and we have set 5 T = since this threshold is typically very 
small. 

In practice the threshold 5 T is small enough to be neglected, and 5\^ have a small impact 
on the b/r ratio, relative to 5^. Therefore, for simplicity, we set <5 r ,i,2 — in this section, and 
we study e {5 ty b,3, tan /3). Furthermore, the dependence of e on 5 t is significantly weaker than on 
5b,3- Hence we use an approximation for 5 t , obtained by taking all superpartners degenerate 
except the Higgsinos, which have mass fi = 500 GeV, and ignoring finite threshold corrections 
to yt- 5t Pprox = 5^ pprox (5s). Below, we will show that deviating from this approximation for 5t 
does not qualitatively change our results. Using this approximation we can draw contours of e 
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in the 5b — £3 plane as shown by the dashed lines in figure |3j with tan/3 = 10(40) in the left 
(right) panel. These dashed lines are model independent; for any values of 5 3i & they yield the 
size of the unified threshold correction necessary for b — r unification. The dashed red line is 
the e = contour that gives precise b — r unification without any unified thresholds, while the 
red shaded area corresponds to unified threshold corrections with magnitude less than 2%. The 
running MS value of the b quark mass is taken to be rrib(mb) = (4.19 + 0.18 — 0.06) GeV at 
3a [29] ; the purple line in figure [3] shows the 2a experimental uncertainty in m&. We determine 
the DR bottom mass, evaluated at Mz, using the 2- loop conversion from MS [3D] arid a 4-loop 
approximation to running a s [31 j . 




Figure 4: The necessary value of tan j3 for t — b — r unification (left), and the GUT-scale Yukawa 
threshold, e (right), as a function of the ?/& and #3 thresholds. At each point, tan/3 (shown as 
green contours to the left) is chosen so that y t = y b at the scale of gauge coupling unification. On 
the right, the GUT threshold, e, is defined by the relation y b (l + e) = y T . As in figure [3j we set 
the gi t 2 and y T thresholds to and St is varied with £3 assuming degenerate scalar superpartners. 

The blue contour in figure [3] shows the relationship between S 3 and Sb that results from 
ignoring all finite thresholds and using the supersymmetric spectrum, described above, with de- 
generate scalars that was used to determine §® pprox (8%) . The superpartner mass scale m increases 
from top to bottom along the blue contour, taking values of (0.1,1,10,100) TeV at the (circle, 
square, triangle, open circle). Note that the finite corrections to 8{ m are negligible since they 
are not tan/3 enhanced, and there are no finite corrections to g$ in the DR scheme [32] • The 
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absence of an intersection between the red and blue lines throughout this range of m and tan /3 
implies that Yukawa unification requires e and/or S[ m to be non-zero. For example, S[ in = 
requires rather large e = 0.2 to 0.3. Alternatively, for a given m and tan/3, one can use this 
figure to determine the required value of S( m to yield e = by reading off the horizontal distance 
between the blue and red contours. As tan/3 decreases, larger values of S( m are requirecQ and 
since S[ m oc tan /3 in this region of tan /3, there is a lower limit on tan /3 that allows e = 0. When 
we consider explicit spectra in section |4j we will find that tan /3 > 10 is necessary when the stops 
and sbottoms are degenerate. 




Figure 5: The relationship between St and 63 (St and S l b ° 9 ) is shown to the left (right) for various 
spectra. For all spectra we fix /i = 500 GeV, assume gaugino unification, and use markers to 
indicate stop masses of mj = 0.1, 1, 10, 100 TeV. For the blue spectrum, the scalars and gluino 
are degenerate with mass rh. For the red (green) spectrum, M 3 = 1/3 (3) rriq. For the purple 
spectrum, M 3 = vn- thl and the first and second generation squarks are 5 times heavier than the 
stop and sbottom. In both plots, the variation among spectra is remarkably mild, justifying our 
simplifying assumption (the blue spectrum) used to relate 5 t and 5 l b ° 9 to £3 in figures [3 and [IJ 

This analysis is extended to the case of t — b — r unification in figure |4j In the left panel, 
for each value of (5b, S3) a value of tan/3 is determined by requiring y t = yb at Mu, yielding the 
green dashed contour lines. In the right panel the dashed lines give the resulting contours for 



2 Larger S[ m is required at lower values of tan /3 because for larger tan ft, is larger and its self-running 
increases yb/yr- 
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e. The e = contour is highlighted in red in both panels, with red shading denoting the region 
with |e| < 0.02. In both panels the blue contour shows the relationship between <5 3 and d>& that 
results from ignoring all finite thresholds and using the same supersymmetric spectrum as in 
figure [3j As rh is increased from 0.1 to 100 TeV, e = requires values of tan/3 from 51 to 48. 
For 5 f b in = 0, tan/3 is near 56. As with lower values of tan/3, sizable values of e and/or 5 b are 
required. 

We conclude this section by testing the robustness of some of the assumptions that we made 
in order to construct figures [3] and |4j In particular, we chose a specific spectrum in order to 
relate 5 t to 5% (which are shown on the right and left vertical axes, respectively) and in order to 
draw the blue curve that relates 5 l b og to 8 3 . The left-side of figure [5] shows S t (5 3 ) for alternative 
spectra having the gluino mass larger or smaller than the degenerate squarks by a factor of 
three, or having the squarks of the first two generations a factor of five heavier than those of the 
third generation. The variations among these spectral lines are remarkably mild, justifying our 
approximate choice of 5 t (5s) in figures [3] and [ij On the right of figure [5J we show 5 t versus 5 l b ° 9 
for light and heavy gluino, and we see that there is no variation in the functional dependence for 
different spectra (the only thing that changes somewhat is the relationship between the scalar 
mass and the position in the S t — 5f, plane). Since 5 l b ° 9 is a fixed function of 5 t , and 5 t is a nearly 
fixed function of £3, these plots together imply that 5 b ° g is a roughly fixed function of £3. This 
means that the position of the blue curve in figures [3] and [4| defined by 5b = S l b ° 9 , also shows 
only mild model-dependence. 

3 Corrections from the Unified Scale 

Over a wide range of parameters, a 10-20% adjustment to yb/y T is required, as shown in figures [3] 
and[4| A key question is how much of this correction originates from the finite contribution to the 
supersymmetric threshold, 5( tn , and how much from the unified scale, e. If the entire correction 
is accounted for by e, then the combination of b — r unification and WIMP LSP dark matter 
does not constrain the scale rh of superpartner masses. 

The unified correction can arise from a tree-level breaking of the SU(5) boundary condition, 
or from loop threshold corrections that depend on the spectrum of the states at the unified scale 

e = e&. c . + (3.1) 

and we consider these in turn. 

3.1 Boundary Condition Corrections 

We focus on unified theories where the dominant contribution to the b and r Yukawa couplings 
are equal, arising from an operator that preserves SU(5). In SU(5) notation, this results from the 
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operator T 3 F 3 H, where T i} F t are the 10 and 5 matter fields for the three generations % = 1,2, 3. 
We allow SU(5) violation in interactions involving (T 2 F 3 , T 3 F 2 , T 2 F 2 ) that lead to the CKM 
element V c b and the mass ratios m s /mb, m ii jm T . From figure [TJ SU(5) violation in these mass 
ratios is required. Given the comparable magnitude of these three quantities - they are all 
within a factor of 2 of 1/30 - we also allow SU(5) breaking in the mixing between the two heavy 
generations. At what level does this SU(5) breaking feed into the b — r boundary condition? 

With the above assumptions, the Yukawa coupling matrices for the heaviest two generations 
in the down (D) and charged lepton (E) sectors of the form 

( B D E 5 C DE 5 \ (0 n \ 

Vo,E = { a , de ! ) V (3.2) 

where B e ^ e ,Cd,e are order unity and 5 = 1/30 governs the scale of (V c b, m s /mb, m^/m T ) . The- 
ories with flavor symmetries acting on T; but not Fi lead to neutrino anarchy, and to hierarchies 
in the up sector that are roughly the square of those in the (D,E) sectors; furthermore, they 
have A' D E of order unity. This leads to large 23 mixing angles on the right-handed quarks 
and left-handed leptons, tan^^ = A' D E , giving a large correction to the boundary condition: 
VblVr = (cos 6 E / cos6d){1 + 0{5 2 )). Hence we assume that any contribution to A' D E larger than 
0{5) is SU(5) invariant 

A' DtE = A + A D>E 5, (3.3) 

where A DjE ~ 0(1). 

After rotating the right-handed states by 6, with tan 9 = A, and redefining parameters, the 
Yukawa matrices take the form 

Vd ' e ={a d>e s i )y- (3 - 4) 

Many theories, especially those based on SO(10), have A = and give this form directly. This 
form is diagonalized by small 0(5) angles giving 

e b .c. = ^(C%-Cl + Al- Al) ~ 10- 2 - 10- 3 . (3.5) 

Hence, in theories of hierarchical charged fermion masses the simplest expectation is €b. c . ~ 5 2 ~ 
Vj ~ 10~ 3 . However, to obtain rris/m^ ~ 3 at the unified scale requires a ratio of Clebsch 
factors of 3 [20J, giving the possibility tb. c . ~ 3 2 <5 2 ~ (m A1 /m T ) 2 ~ 10~ 2 ; hence the range given 
in equation (3.5). Such a Clebsch might originate in the 22 entry of the Yukawa matrix and 
be transferred to C EiE via a large 23 rotation due to A ~ 0(1), or it might be present in the 
original 23 entry. 

A wide class of unified theories, having Yukawa matrices of (3.4), lead to e^.c. < 0.01. In these 
theories b — r unification requires (10-20)% threshold corrections from either supersymmetric or 
unified scales. 
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It is important to note that theories having an SU(5) breaking contribution to the 33 entry of 
Ud,e °f order 6 will typically give e&. c . ~ 5 ~ 3%. If this is further enhanced it could allow b — r 
unification without threshold corrections [33J. Consider an SO(10) theory where the Yukawa 
interactions involving the third generation are 

2/^3^3 + y'(cos9i/; 3 + smOfo) 2 </>', (3.6) 

where the 16-plet spinors ipi contain Tj and Fj. The multiplets and 0' are 10 and 126 dimen- 
sional and contain components of the MSSM Higgs doublets D x a H a , 0' D x' a H a , a = u,d. A 
weak triplet in 0' acquires a vev leading to type II seesaw neutrino masses, with 9 ~ 45° the 
atmospheric neutrino mixing angle. The coupling y 1 must be chosen such that the weak doublet 
vev in 0' yields V^, implying that the , 03'030 / interaction leads to 

e b ,. =4V cb -^~0.1-^- (3.7) 
1 — x 1 — x 

where x = x u x' d /x^x' u . 

3.2 Threshold Corrections 

Contributions to the unified threshold correction e t h are typically comparable to those that 
correct the gauge couplings, e a , since they both arise at 1-loop, involve the same couplings 
and the same logarithms of mass ratios. The success of precision gauge coupling unification 
then implies that |e^| < 0.02, so that it can be neglected. However, there are exceptions to 
this situation and we consider two in this sub-section. Some couplings contribute at 1-loop 
to logarithmic Yukawa thresholds but not gauge thresholds. The top Yukawa coupling is an 
example, but is usually not large enough to give significant contributions to e t h- However, if 
tan/3 is sufficiently reduced these corrections can become important. Secondly, in orbifold grand 
unified theories the KK towers at the unified scale do make substantial corrections to gauge 
coupling unification, so it is important to understand how they may affect Yukawa unification. 

3.2.1 From the Top Coupling 

At low tan /3 the top Yukawa coupling increases, changing the evolution of the b quark Yukawa 
coupling in the direction of improving b — r unification. This can allow the squark mass scale 
to increase well beyond 10 TeV while keeping Higgsinos at the TeV scale for dark matter. An 
example of this situation is shown in the left panel of figure [6] where it is compared to the case 
of precision Yukawa unification at high tan /3. The blue curves show the evolution of y^ with 
tan/3 = 20 for 5 f b in = (dashed) and 0.23 (solid). The latter case leads to precision unification 
with the t coupling, shown by the red solid curve. The squark masses are taken to be 2 TeV 
and this is an illustration of the precision Yukawa unification discussed throughout this paper. 
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Figure 6: The left panel shows the evolution of the b Yukawa coupling in two schemes. The 
blue curves are for tan (3 = 20 and squarks at 2 TeV, without (dashed) and with (solid) finite 
supersymmetric threshold corrections. The green dotted curve has a split spectrum with squarks 
at 10 3 TeV and tan/3 1. The red solid curve shows the evolution of the r Yukawa coupling 
in both schemes. In each scheme the curves are normalized by the value of the unified Yukawa 
coupling. In the right panel the solid blue curve shows the required e for Yukawa unification in 
the scheme with the split spectrum as tan (3 is varied, and the dashed purple line shows that the 
unified top Yukawa is rapidly increasing at low tan /3. 

The dotted green line illustrates the evolution of the b Yukawa coupling for an entirely 
different possibility: split supersymmetry with squarks at 10 3 TeV. The fermionic superpartners 
are at the TeV scale with /x = 0.5 TeV, M 3 = 2 TeV and gaugino mass unification. Since the A 
parameter breaks R symmetry it is expected to be of order the gaugino masses and makes no 
relevant radiative contribution to the Higgs mass. Indeed, we take tan (3 ~ 1 so that the entire 
Higgs mass arises from RG scaling of the Higgs quartic between the top squark mass and the weak 
scale, requiring a stop quark mass of about 10 3 TeV for a 125 GeV Higgs mass. In this example 
all squark masses are set at 10 3 TeV. With such heavy squarks there is an important logarithmic 
supersymmetric threshold to b/r, which works against Yukawa unification [25]. Hence the need 
for very low tan (3 leading to a large y t . 

The effect of increasing y t is seen by comparing the dashed blue and green curves. The 
difference in the running of these curves is most pronounced near the unification scale, where 
the top Yukawa coupling in the tan (3 = 1 case is getting very large. The dashed purple curve in 
the right panel of figure [6] shows y 2 /16ir 2 at the unification scale; it is less than 0.01 at large tan (3 
and rises very rapidly to over 0.2 by tan /3 = 1.03. The corresponding reduction in the size of the 
required correction e at the unified scale is shown by the solid blue curve. Indeed at tan/3 = 1.03 
no correction is needed. However, this cannot be viewed as precision Yukawa unification since 
unknown threshold corrections of order y 2 /16n 2 = 0.2 are expected at this point, even if they are 
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not logarithmically enhanced. Perhaps the best one can do is near tan/3 = 1.05, where e = 0.06 
is required and corrections of order rfi/lQir 2 = 0.06 are expected. More generally the situation 
is clear: Yukawa unification at low tan does not have a precision that is controlled by the low 
energy theory, rather the required 10-20% correction is coming from the unified scale where the 
top Yukawa coupling is approaching its Landau pole. 

3.2.2 From KK Towers in Orbifold GUTs 

Precision gauge coupling unification can occur in theories where the unified symmetry is realized 
in higher dimensions but is broken by the compactification to 4d. The higher dimensional 
symmetry forbids Yukawa couplings from appearing in the bulk, and hence Yukawa unification 
can occur in such theories if the operator generating the b and r Yukawa couplings is located on 
a fixed point that respects SU(5). Such Yukawa couplings do not have threshold corrections from 
KK modes because of locality. However, gauge couplings do receive corrections from KK towers, 
and this can alter the scale of gauge coupling, and therefore Yukawa coupling, unification. 

In the simplest 5d orbifold GUT [34J KK corrections improve gauge coupling unification 
compared to 4d unification, and reduce the unification scale to 10 15 GeV. The b and r couplings 
must therefore unify at this scale with e = 0, as shown for tan /3 near 50 by the green dashed line 
in figure [7J Similar figures apply for lower tan/3. Thus b{ m ~ 0.1 — 0.15 is an absolute necessity 
in such theories. 

3.2.3 Summary 

We have argued that there is a wide class of unified theories having small boundary condition 
and threshold corrections to the bjr Yukawa ratio at the unified scale, |e| < 0.02. Our analysis 
of supersymmetric threshold corrections in section [4] and the resulting collider, dark matter and 
B meson physics of section [5] applies to these theories. Nevertheless, it is possible to construct 
theories where e is much larger. Such unified theories are required if supersymmetry is broken 
to the SM at the unified scale, when e ~ 0.6 is needed. 

4 MSSM Spectrum from Yukawa Unification and the 
Higgs mass 

In this section, we determine which supersymmetric spectra allow for third generation Yukawa 
unification and a Higgs mass of 125 GeV. In order to remove dependency on unknown GUT 
physics, we simply assume that the GUT-scale threshold corrections to yb/y T are small enough 
to be neglected, |e| < 0.02. This assumption, which is motivated by the observation that the 
gauge thresholds are also this small, applies to a broad class of GUT models as described in 
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Figure 7: Yukawa unification in an orbifold GUT with compactification scale M c = 10 15 GeV. 
The black contours show the b — r unification quality, e, evaluated at the compactification scale, 
as a function of the y^ and g 3 threshold corrections. At each point, tan/3 is chosen so that y t = 
at the compactification scale, and for simplicity we set the gi^ and y T thresholds to 0. The green 
shaded region denotes where |e| < 0.02 in the orbifold GUT. For comparison, the red shaded 
region shows where |e| < 0.02 in the 4D GUT, as in figure [ij 



section [3j Then, the success or failure of Yukawa unification depends entirely on the weak-scale 
SUSY spectrum, which must be chosen so that S( m ~ 0.1 — 0.2, as we saw in section [2j The 
program of this section is to determine which spectra satisfy this criterion. We will also impose 
the requirement that the LSP abundance not overdose the universe, which when taken together 
with Yukawa unification implies, as we will see, an upper bound on the scalar masses. 

The most important threshold corrections for yb/y T are the corrections to y^, yt, and g 3 , as we 
discussed in section M For this section, we use the complete 1-loop SUSY threshold corrections, 
as computed by Ref. [32]. The qualitative behavior of these thresholds is simple to understand. 
The #3 threshold is defined to be logarithmic only, in the DR scheme, and the yt threshold is 
dominated by the log contribution. Therefore, both St and S3 are negative, with size determined 
logarithmically by the overall SUSY mass scale. The yj, threshold is a sum of the negative 
logarithmic piece and a tan/3 enhanced finite piece, S[ in . As we found in section [2J the b and 
r Yukawas unify when S( m ~ 0.1 — 0.2, where S[ m is dominated by the gluino and chargino 
diagrams shown in figure |2j which is approximately given by equation 1.1 The most important 
SUSY parameters for determining S[ m , and therefore whether or not the Yukawas unify, are, 



16 



mg 3 , mu 3 , m D3J A t , M 3 , /i, tan/?. 



(4.1) 



Note that throughout this paper, we use the SUSY Les Houches Accord (SLHA) conventions for 
the signs of these parameters [35, 36J. 

We choose to take seriously the hint, from ATLAS pQ and CMS [2], that m^ w 125 GeV. 
We will check to see that Yukawa unification can be made consistent with a Higgs boson of this 
mass. Conveniently, the requirement that the correct Higgs mass is generated radiatively by 
the SUSY spectrum allows us to reduce the SUSY parameter space by one dimension. Recall 
that the tree-level Higgs mass is given, in the MSSM, by m\ = cos 2 2/3 m\ ~ m 2 z (where we 
have taken the decoupling limit, 3> mz, and used the fact that in the parameter space of 
interest for Yukawa unification, tan/3 is large enough that cos 2 2/3 « 1). Radiative corrections 
are required to raise the Higgs mass from mz to 125 GeV. The leading 1-loop correction is given 
by, 



cos 2 2/3 m% + 



4tt 2 



4 2 

Vt v 



log 



+ 



X 



mr 



in 



1 



Xf 



12m? 



125 GeV, 



(4.2) 



where v 



17 '4 GeV, m\ 



mQ 3 mu 3 , and X t = A t — /// tan j3 determines the amount of left-right 



stop mixing. By solving equation 4.2, we can determine A t for a given value of mi 



In practice, the RG improvement of equation 4^2 must be taken into account to achieve 
accurate results. In order to determine the Higgs mass as a function of the soft parameters, 
we use a version of Suspect designed to correctly handle heavy scalars, provided to us by the 
authors of Ref. |37J. In this code, the Higgs quartic coupling A is run from the weak scale to the 
stop mass, mj, using 1-loop RG and 1-loop thresholds computed both at the weak scale and the 
stop mass scale. The result for tan/3 = 50 is shown in figure |8j as a function of mi and X t /mi. 
The red band shows the uncertainty coming from the top mass, mt = 173.2 ± 0.9 GeV (38]. For 
small values of the A-term, X t ~ 0, we require mi ^4 — 5 TeV. The smallest value of the stop 
mass corresponds to the regime known as maximal mixing, where X t is chosen to maximize the 
second term in the brackets of equation 4.2, X t /mi = ±\/6. Here, we find m t ~ rs 900 GeVj^] 



Moving to larger values of \X t \, the scalars drift up to heavier masses. For X t /mi = ±vL2, the 



X t dependent terms of equation 4.2 cancel against each other. For larger values of the X t 



dependent terms are negative, and can be chosen to cancel against the log, allowing for heavier 
scalars. 



3 Note that near maximal mixing, the 2-loop contributions to mh are important for determining the stop mass 
and there remains a large amount of uncertainty on the necessary stop mass, with different estimates ranging 
from mj ~ 700— 1000 GeV, see Ref. [3]. We do not include the 2-loop effects in this paper, since the public codes 
implementing the 2-loop Higgs mass SO] are unreliable in the limit of heavy scalars. 
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Figure 8: The values of the stop mass, m~ = mQ 3 mu 3 and stop mixing, X t , that lead to 
rrih = 125 GeV, with tan (3 = 50. The red band shows the uncertainty coming from the top 
mass, m t = 173.2 ± 0.9 GeV. Here the Higgs mass is calculated at 1-loop, using a numerical 
code [37] that allows for heavy stops by re-summing \og(mi/mz)- 

We have seen that rrih « 125 GeV potentially allows for very heavy scalars for special, large, 
values of A t . In this part of parameter space, the large A t leads to the presence of a lower 
energy charge and color breaking vacuum [41], and our vacuum is metastable. The criterion for 
global stability is A\ + 3/i 2 < 3(mn 3 + m^ 3 ) [42]. If we take rriQ 3 = mu 3 = and assume that 
\i <C At, mf, then metastability corresponds precisely to X t past maximal mixing, \X t \ > v^6 m t ~, 
where we have used that X t ps A t at large tan/3. Splitting m,Q 3 from mu 3 , while keeping the 
Higgs mass fixed, has the effect of increasing the globally stable region. For X t past maximal 
mixing, one must check whether or not the lifetime of our vacuum is long enough to support 
the present Age of the Universe. The numerical analysis of Ref. [43] found that our vacuum is 
sufficiently long-lived for Af + 3fi 2 < 7.5(mg 3 + m^). For degenerate mg 3 and mjj 3) this would 
correspond to X t /m t - < -\/l5, or mi < 54 TeV for rrih = 125 GeV. Unfortunately, the numerical 
analysis of Ref. [43] was only performed for < 2 TeV and needs to be extended to include 
the part of parameter space relevant for rrih ~ 125 GeV. We leave this for future work, and for 
now we tentatively consider the entire parameter space that gives rrih = 125 GeV. 

We show our primary result in figure [9] the SUSY parameter space that allows for precise 
b/r unification. As discussed above, we have related and A t by imposing rrih ~ 125 GeV. 
For simplicity, we take degenerate stops and sbottoms (we will relax this assumption below), 
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Figure 9: The SUSY parameter space with third generation Yukawa unification, assuming 
degenerate stops and sbottoms rriQ 3 = m Ua = m D3 , for tan/3 = 50,30, 15. If tan/3 = 50, then 
t — b — r unification is possible; for the lower values of tan /3 only b — r unification is possible. We 
determine the stop mass as a function of the stop mixing, X t , by demanding = 125 GeV, as 
shown in the lower row of plots. In the upper plots, we show the value of /i that gives Yukawa 
unification as a function of X t and M 3 . The green shaded area has // < 1 TeV, allowing for a 
bino-Higgsino LSP whose energy density does not overdose the universe. If the LSP is wino-like 
and constitutes all of DM, then M2 ~ 2.7 TeV and we require \i > M2, which is the area outside 
the purple contour. 

m| = rriQ 3 = m\j 3 = mj^ 3 . The top row of plots show the value of fi that leads to b — r 
unification, as a function of X t jm- t and M 3 , for tan/3 = 50,30, 15. The first value, tan/3 = 50, 
leads to full t — b — r unification, as we saw in figure [4j while the lower values of tan /3 allow 
for b — t unification only. The bottom row of plots shows the corresponding value of m 4 -, as 
a function of the X t /mi axis, that gives mn = 125 GeV. The remaining SUSY parameters 
have sub-leading effects on Yukawa unification, but to be definite, we fix Mi j2 using the gaugino 
unification relation, 6M\ w 3M2 ~ M3, all of the other sfermions are fixed to mf, and we take 
rriA = 3 TeV. 

The interplay between Yukawa unification and the cosmological LSP abundance is key. As 
stated in the introduction, we consider high-scale mediation and conserved R parity so that the 



19 



LSP is stable and has an abundance in standard cosmology given by thermal freezeout. The 
shaded green region with \i < TeV is the one of primary concern to us, since it corresponds to 
the simplest DM possibility of bino/Higgsino LSP. Note that a pure Higgsino has the correct 
abundance, Q c d m h 2 = 0.11 [29J, when fi « 1 TeV. A bino-Higgsino admixture does not overdose, 
fijY < Qcdm, as long as /i < 1 TeV and \i < Mi, and a bino-Higgsino admixture typically has 
the correct abundance to constitute all of dark matter in the so-called well-tempered regime, 
fi ~ Mi [27]. Since the precise value of M\ has little effect on Yukawa unification, anywhere 
within the green region can be made to avoid overclosure simply by choosing M\ > fi. 

There are several other possibilities for SUSY dark matter. A bino/Higgsino admixture can 
have fi > 1 TeV and still avoid overclosure if DM co-annihilates with a slepton or squark or if 
DM annihilates through a heavy Higgs pole. Both of these options necessitate extra light scalars 
and require a tuning among masses to avoid overclosure. We will not consider co-annihilation 
or heavy Higgs pole further in this paper. Another option, that does not necessitate tuning 
parameters to avoid overclosure, is that the LSP is wino-like, with M2 < Mi, /i. Then the 
correct abundance for the wino to be all of DM is achieved for M 2 ~ 2.7 TeV (SJ US], an d 
lighter winos would contribute a subdominant portion of the energy density. If we assume that a 
wino LSP contributes all of DM and demand Yukawa unification then we must have \x > 2.7 TeV: 
this boundary is shown SIS db purple line in figure [9] 

Returning to the simplest possibility: that DM is a bino/Higgsino admixture, we see from 
figure [9j that avoiding overclosure, \i < 1 TeV, leads to an upper bound on the stop mass when 
we demand Yukawa unification. For tan = 50, as is appropriate for t — b — r unification, we 
find that m t ~ < 2.8 TeV is required when A t and M3 have opposite signs, which corresponds to 
the right half of the plane shown in the figure. Here, the gluino and chargino diagrams have 
opposite signs: see equation |1.1[ We note that a relative sign between A t and M 3 is favored 
by the renormalization group because M 3 contributes radiatively to A t , driving it to have an 
opposite sign, 

^-^(^3 + 6|rfA + W^ + ...). (4.3) 

In order for A t and M3 to have the same sign, there must be a large A-term in the UV to 
compensate the radiative contribution from M3. In this case, the two diagrams add and heavier 
stops are possible, mi < 8500 GeV. Since the thresholds are proportional to tan/3, lowering 
tan/3 has the effect of shrinking the allowed parameter space and necessitating even lighter 
stops (also, recall from figure [3] that a larger 5( m is required at lower tan/3). For tan/3 = 30, 
we require mi < 1400 GeV when A t and M 3 have opposite signs, and m 4 ~ < 4600 GeV when 
they have the same sign. For tan /3 = 15, the region where A t and M3 have opposite signs is 
removed completely, and we require mi < 2000 GeV when the signs are the same. We find no 
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Figure 10: The allowed parameter space for Yukawa unification, as in figure [9j except relaxing 
the assumption of degenerate stops and sbottoms. On the left, we consider lighter right-handed 
sbottoms, with mo 3 = 1/3 and m 4 ~ = mg 3 = mjj 3 . On the right, we consider heavier right- 
handed sbottoms, m£> 3 = 3m ( -. Lighter sbottoms enhance the size of the gluino diagram of 
figure [2j leading to a larger region with Yukawa unification, while heavier sbottoms deplete the 
gluino diagram and shrink the allowed region. 



b — t unification at all when tan f3 < 10. We note that with our sign conventions, nearly the 
entire parameter space with Yukawa unification and ji < 1 TeV is where M 3 is negative, so we 
have only shown this sign choice in figure [9] 

So far, we have assumed that the stops and sbottoms are degenerate. We relax this assump- 
tion in figure 10 On the left, we show a case with lighter right-handed sbottom, m£, 3 = 1/3 mi, 
where mj = mg 3 = mu 3 . On the right, we show the case with heavier right-handed sbottom, 
m£, 3 = 3m ( -. The effect of splitting the stops from the sbottoms is to change the relative size of 
the gluino diagram, which has sbottoms propagating in the loop, to the chargino diagram, which 
has stops propagating in the loop. When the sbottoms are lighter, as in the left diagram, the size 
of the gluino diagram is enhanced, leading to a larger region with Yukawa unification. In this 
case, stops as heavy as 14 TeV allow for Yukawa unification when A t and M 3 have the same sign. 
However, this implies that the sbottom is lighter than 5 TeV, such that the requirement of one 
(somewhat) light scalar has not been circumvented. When the sbottoms are heavier, as in the 
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right diagram, the gluino diagram is depleted. In this case, we find no Yukawa unification when 
A t and M 3 have opposite sign, and when they have the same sign we require < 5200 GeV. 



5 Flavor and Dark Matter Phenomenology 

In this section we discuss possible experimental signatures of supersymmetry with Yukawa unifi- 
cation. We found in section [1] that most of the parameter space of Yukawa unification (plus dark 
matter) has stops and sbottoms lighter than several TeV and therefore potentially accessible to 



experiments. We discuss the flavor signals and constraints in section 5.1 As we will see, the large 



values of tan /3 necessary for Yukawa unification lead to observable signals in I?-meson decays. 



In section 5.2, we discuss the phenomenology of dark matter within our parameter space. 

Before we move on to discus flavor and dark matter, we briefly consider the exciting prospect 
of discovering superpartners at the LHC. Yukawa unification implies that the stop and sbottom 
must be lighter than several TeV and therefore may be produced at the LHC. The observability 
is highly dependent on the precise spectrum. The masses of the superpartners will determine 
the production cross-section, and their spectrum will determine which decays dominate and 
therefore whether or not low-background final states (preferably with leptons and fe-jets) are 
populated. 



As a simple example of the overall SUSY production rate, figure [TT] shows the cross-section 
as a function of the gluino mass and a common squark mass, at center of mass energies of 13 and 
14 TeV. If a search in a low-background channel has an efficiency of ~ 10%, then 50 events before 
cuts are more than enough for a discovery. For degenerate gluino and squarks, this corresponds 
to a 100 fb _1 reach of 2.4 (2.3) TeV at 14 (13) TeV. It is important to keep in mind that discovery 
is not guaranteed because while Yukawa unification requires a light stop and sbottom, it does not 
require a light gluino (see section [4j and the masses of the first and second generation squarks 
are unconstrained. Still, assuming a simple scenario for supersymmetry breaking, the stop and 
sbottom masses are probably near the masses of the other colored sparticles. Therefore, we find 
it highly encouraging that the eventual LHC energy will be enough to probe most superpartner 
scales relevant for Yukawa unification. 



5.1 5-meson decays 

In this section we consider the flavor violating signals of SUSY with Yukawa unification. We 
found in section [4] that the favored parameter space of Yukawa unification has stops and sbottoms 
with masses in the ~ 1 — 10 TeV range. For superpartners at these scales, generic squark soft 
masses are forbidden by AF = 2 processes such as K — K mixing. Therefore the soft masses 
must possess a special flavor structure, such as Minimal Flavor Violation (MFV) [4*81 1491 150] or 
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Figure 11: The total colored sparticle production cross-section, at LO, as a function of the 
gluino mass and a common squark mass [4tj| . The red area is already excluded by ATLAS, 
assuming the squarks and gluinos decay directly to a massless neutralino [27] . 

a £7(2) flavor symmetry [3Tj 132"] 153"] (or U(2) 3 [3D). Even if we do assume MFV, there are B- 
meson decays that receive tan (3 enhanced contributions [33] . This is highly relevant for Yukawa 
unification, because we found above that t — b — r unification requires tan /3 ~ 50 and b — r 
unification plus dark matter requires tan j3 > 10. For now, we assume that CP phases are small 
and only consider flavor violating effects. We will comment on the limits on CP violation at the 
end of this section. 

The most constraining processes are b — > s^y and B s — > / u + . The leading SUSY diagrams, 



at large tan/3, arise from gluino and Higgsino exchange, and are shown in figure 12 There is 
a tight relationship between the amplitudes for these processes and the threshold correction to 
the bottom mass, which is generated by similar 1-loop diagrams, as can be seen by comparing 



figures [2] and 12 Limits from 5-meson decays are discussed in many papers on Yukawa unifica- 



tion, see for example Refs. [21] |23l E3J [56]. These references all consider unified scalar masses at 
the GUT scale. Our approach differs because we allow for general soft terms at the weak scale, 
without specifying the high-scale boundary condition. 

We now consider the flavor structure of these diagrams. We work in the basis where the 
fermion masses and squark soft masses have been diagonalized, with conventions described in 
appendix[B] In the limit of small left-right mixing, this corresponds to the squark mass eigenstate 
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Figure 12: The leading, tan/3 enhanced, diagrams contributing to b — > 57 and B s — >■ 
are shown above and below, respectively. There is a close relationship between the size of these 
diagrams and the finite threshold correction to the bottom mass (see figure [2]). In the SUSY-CKM 
basis, the leading contributions of the gluino exchange diagrams on the left are proportional to 
left-left down squark mass mixing, (5d L ) 32 , while the Hig gsino exchange diagrams on the right 
are proportional to V ts , coming from the Higgsino vertex. 

basis, and all flavor violation can be written in terms of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix and the gluino vertices, (W q ) v ' qjqjg (see Ref. [57] and appendix Sj). Small left-right 
mixing is a good approximation in most of our parameter space, because the squarks are heavier 
than 1 TeV and left-right mixing is suppressed by v/rh. 

First consider the Higgsino diagrams contributing to b — > 37 and B s — > // + /i~ . The diagrams 
for the two processes are proportional to the same flavor factor, but involve different loop func- 
tions. Assuming that the Higgsino diagrams are dominated by stop exchange, because of the 
A-term insertion, 

iM k oc y b y t A t V ts (W£W^)* « y b y t A t V ts , (5.1) 

where the second step assumes small mixing, so that the 3-3 entries of W u ° and W d are close 
to 1. We see that the Higgsino diagrams have an irreducible contribution proportional to Vt s , 
arising from the Higgsino vertex connecting a strange quark and top squark. 

Next we turn to the gluino diagrams. The leading contribution has the flavor structure, 

iM s cx (w dT 2j P $ W d * 3 ) y b (W dc l P 3 C W?;) » y b P< (w d 2 P 3 + W%P 2 ) (5.2) 
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where P, (P£) denotes the propagator of (<i 3 ) and the second step assumes 1-3 mixing is small 
enough to be neglected. Note that this diagram vanishes for degenerate squarks because of a 
super-GIM mechanism that follows from the unitarity of W d . For non-degenerate squarks, the 
diagram is proportional to W^, which is an unknown factor of order Vu in theories with MFV. 

It is also common to describe the flavor violation in the so-called SUSY-CKM basis, where 
the entire quark superfields are rotated into the basis where the fermion masses are diagonal. In 
this basis, it is useful to work in the mass-insertion-approximation when flavor violation is small, 
in which case the leading contribution from the gluino diagram is proportional to the left-handed 
2-3 mixing, (5% L ) 32 . The factor (5% L ) 32 is related to W 23 by equation B.6, and is expected to be of 
order V ts in theories with MFV. In general, the gluino diagrams also contain terms proportional 
to the right-handed squark mixing, {S RR ) 32 . This mixing is suppressed in theories with MFV 
and here we assume that it is subdominant to the contribution from left-handed squark mixing. 

We now describe the experimental constraints, and reach, in some detail. We will first discuss 
the constraint coming from b — > S7, which has already been measured and agrees with the SM 
prediction. Then we will discuss B s — > which has not yet been measured and, as we will 

see, may have a rate that differs significantly from the SM prediction. 

The branching ratio for b — > 57 has been precisely measured by the B-factories to be (3.55 ± 
0.26) x 10~ 4 [58J , which agrees with the NNLO SM prediction of (3.15 ± 0.23) x 10" 4 [59]. This 
presents a serious constraint on theories with large tan/5, even for squarks in the several TeV 
range, because the amplitude of the SUSY contribution is proportional to tan (3. Figure 13 shows 
the allowed region for Yukawa unification for tan /3 = 50, as in figure [9] except now shown as a 
function of the stop mass and gluino mass. As before, the blue contours indicate the values of 
\i necessary for Yukawa unification, and the shaded green region has fi < 1 TeV and is therefore 
compatible with a bino/Higgsino LSP that does not overdose the Universe. The red shaded area 
would be excluded at 2a by b — > 57 if we were only to consider the irreducible contribution from 
Higgsino exchange, by setting (5% L ) 32 = 0, which turns off the gluino diagram. Note that for our 
numerical results, we for simplicity use the mass-insertion-approximation results of Ref. [60] . 

If the Higgsino diagram were the end of the story, we would conclude that a large portion of 
the Yukawa-unified parameter space is already excluded. However, it is important to include the 
gluino contribution, which depends on the unknown coefficient, {5% L ) 32 . We show the branching 



ratio of b — > 57 as a function of (^ L ) 32 i n figure 14 for the points indicated by black, blue, 



and purple dots in figure 13 The shaded area of the plot is excluded at 2a and we see that 
each point is safe for a wide range of values of {S^ L ) 32 . We find that no fine-tuning is necessary 
because the Higgsino and gluino diagrams are of similar size and their sum is easily consistent 



with the constraint. In figure 13, the red contours indicate the minimum absolute value of 
(5^ L ) 32 necessary for b — > 57 to be consistent with observation, in units of | V* s | . In most of the 
parameter space, a value of 0(y ts ) is sufficient. 
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Figure 13: The allowed region for Yukawa unification with tan/3 = 50, with the limit from 
b — > S7 overlaid. The y-axis is the gluino mass and the x-axis corresponds both to the stop 
mixing, X t and the stop mass, which is non-monotonic in X t because after fixing = 125 GeV, 
the stop mass increases moving away from maximal mixing, X t = ±\/&mi (see figure ^J. As 
in figure |9j the blue contours indicate the value of /i necessary for Yukawa unification and the 
shaded green region is consistent with bino/Higgsino dark matter, /i < 1 TeV. The red region 
is excluded, at 2a, by b — >• 37 if 2-3 down squark mixing is ignored, (^ L ) 32 = 0, such that 
only Higgsino exchange contributes (we take uia = 10 TeV so that charged Higgs exchange is 
decoupled). The red contours show the minimum value of |(^ L ) 32 |, in units of \V ts \, necessary 
to bring b — > 57 in accord with the observed branching ratio. The black, blue, and purple points 



are considered in more detail in figures 14 and 15 



We make a few technical comments about b — > 57 before moving on. In addition to the gluino 
and Higgsino exchange diagrams, there is a also a diagram with wino exchange, where the wino 
mixes with a Higgsino. This diagram is proportional to the same flavor structure appearing in the 



first factor of parentheses in equation |5.2| The wino exchange is typically subdominant to gluino 
and Higgsino exchange, but we have included it in our numerical analysis by setting the ratio 
M 3 : M 2 = 3, as is approximately satisfied by spectra with gaugino unification. Our answers are 
not very sensitive to the choice of M 2 . A second comment is that we have been focusing on the 
allowed region of b — > 57 where the SUSY contribution is small. There also used to be an allowed 
region where the SUSY contribution has the opposite sign of the SM contribution, but where 
their sum happens to have the same magnitude as the SM contribution. This finely-tuned option 
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Figure 14: The branching ratio for 6—7-57 normalized to the SM prediction of (3.15 ± 0.23) x 
10~ 4 [59], as a function of (S^ L ) 32 in units of \V ts \. The left (right) side of each plot corresponds 
to {S^ L ) 32 < (> 0). The black, blue, and purple curves correspond to the corresponding points 



in parameter space denoted on figure [13} The shaded gray region is excluded by more than 
2a. Notice 
wide range 
decoupled. 



2a. Notice that for each point, the branching ratio of b — > 57 is consistent with the limit for a 

)32- 



wide range of values of (8% L ) 32 . We have taken = 10 TeV so that charged Higgs exchange is 



is now excluded by LHCb measurements of b — > X s l + l~ [61]. In the region we are interested in, 
the SUSY contribution to b — > 57 is small, and in this limit b — > X s l + l~ is not relevant. 

More interesting than b — > S7, which is already observed to agree with SM, is the potential 
to see an observable deviation in the rate for B s — > ji + /i~ . This process is predicted to have the 
branching ratio (3.2 ± 0.2) x 10~ 9 in the SM [62]. The experimental limits have been rapidly 
approaching this target, and LHCb now sets a 95% limit of 4.5 x 10~ 9 [63]. CMS is close behind 
with a 95% limit of 7.7 x 10~ 9 [61]. It is likely that we will know very soon whether or not 
this rate agrees with the SM prediction or shows an observable enhancement or, in the case of 
destructive interference with SM diagrams, depletion. The decay is generated by the 4-fermion 
operator, 

C eff D C s m b s L b R fifx, (5.3) 

and also the operator with a 75 between the muons, and the corresponding operators derived by 
flipping parity, L <H- R. 

The leading SUSY contribution is mediated by a heavy Higgs and has an amplitude pro- 
portional to tan 3 /3 |65J. This contribution follows directly from the flavor- violating analogue of 
the finite threshold correction to the bottom mass, which induces, in mass eigenstate, an sib^H 
coupling. Therefore, the SUSY contribution to the coefficient of this operator can be written 
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Figure 15: The reach for discovering a modification to B s — > fi + {i , for the two points in 



parameter space denoted by black dots on figure 13, as a function of the heavy pseudoscalar 
Higgs mass, m^, and (S^ L ) 32 i n units of \V ts \. The shaded gray region is excluded, at 2a, by 
b — > 57. The purple contours indicate the branching ratio of B s — > fi + fx~ in parts per billion. 
The red region is excluded at 2a by LHCb with 1 fb _1 , Br^^-n^- < 4.5 x 1CT 9 . The blue (green) 
regions have a branching ratio that is smaller (larger) than the predicted SM cross-section of 
(3.2 ± 0.2) x 10~ 9 [62J by 3a in the region of light shading and 5a in the region with darker 
shading. 



directly in terms of the finite bottom threshold, dressed by the appropriate flavor factors, 

\^ L ) 32 si + v;X], (5.4) 



(jSUSY _ m H tan2 $ 



'S ~ 9 2 



where 5? and 5^ denote the gluino and Higgsino contributions to the finite bottom threshold, 



which are given by the two terms in equation 1.1 The size of the SUSY contribution is highly 



sensitive to tan /3 and to the value of the heavy Higgs mass. 



The current limit, and reach, for discovering B s — > are displayed in figure 15 



as a 



function of the squark mixing, (S% L ) 32 , and the pseudoscalar Higgs mass, m^, for the two points 



denoted by black dots in figure 13, The shaded gray regions are excluded by b — > 57 (note that 
we have included the charged Higgs contribution, which effects the boundary of the excluded 
region for rriA < 1 TeV). The shaded red area is excluded by the LHCb limit on Br Bs _ 5 . At + /J -. 
In nearly the entire parameter space this limit is significantly stronger than the direct LHC 
limit on heavy Higgses. The strongest direct limit is set by CMS on H — > t + t~, requiring 
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Figure 16: The limit on b — > 37 and the reach for B s — > fi + fj, are shown for tan (3 = 20 on the 
left, and right, respectively. The left plot corresponds to the region with X t < and all contours 



are as in figure 14 The right plot corresponds to the point shown as a black dot in the left plot, 



and all contours and shading are as in figure [15 



rriA > 480 GeV when tan (5 = 50 [66J . The (dark) blue region denotes values of the branching 
ratio of B s — > fi + {i~ 3a (5a) lower than the SM prediction]^] We see that a large range of rriA 
leads to an observable depletion for the point with X t < 0. The green regions denote values 
that are enhanced by 3a and 5a. This reflects a much smaller portion of the available parameter 
space because of the proximity of the stringent LHCb limit to the SM prediction. 

So far we have assumed tan /3 = 50 in this section, as is appropriate for t — b — r unification. 
But we found that b — t unification is possible as long as tan f3 > 10, and one may wonder 
how lowering tan /3 impacts the limits from 5-meson decays. In order to explore this, we choose 
the example value of tan/3 = 20, and in figure [16] we show the b — > S7 limit on the left, for 
X t < 0. Interestingly, we find that the b — > S7 limit is stronger at lower values of tan 0: a 
larger fraction of the region with fi < 1 TeV is excluded by the Higgsino-mediated diagram and 
therefore requires some cancellation between the Higgsino and gluino diagrams. This might be 



Note that the 3cr and ha regions of figure 15 utilize the theory uncertainty of Ref. [62], where the Fb b 
uncertainty has been removed by assuming that the measurement of AAIg s is SM-like. Any deviation from the 
SM prediction needs to be treated with care because new physics can contribute to AMb s - In SUSY theories 
with MFV, new contributions to AMb s can be neglected, relative to new contributions to B s — > because 
they rely on R-R squark mixing, which is small. 
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seem surprising since the amplitude for b — > S7 is proportional to tan f3. However, so is the finite 
threshold to the bottom mass, and the b — > sj amplitude scales roughly as 5( m /m^. For a 
given value of /1, as tan (3 is reduced the stop and sbottom masses must be lowered to keep S[ in 
large enough for b — r unification, resulting in a stronger b — > 57 limit. In addition, as we lower 
tan/3, a slightly larger S( m is required (see figure p|. Still, we find that only a mild cancellation 
is required between the Higgsino and gluino diagrams to bring b — > 57 in accord with the limit, 
for most of parameter space, even as tan f3 is lowered. 



We now turn to the B s — > jj, reach at lower tan j3, which is shown to the right of figure 16 



for tan (3 = 20 and the point in parameter space denoted by a black dot on the left side of the 
figure. We see that the LHCb limit, and region with a reach to observe an enhanced or depleted 
branching ratio, are both pushed to lower values of relative to the tan = 50 scenario. This 



is because the amplitude for B s — > /i + scales as tan 2 (3 x S[ m /m\ (see equation 5.4), and the 
lower value of tan f3 is compensated by a lower value for m^. Note that a large fraction of the 
allowed values of vtia still lead to an observable effect since the present limit is also reduceoQ 

We now discuss the relationship between Yukawa Unification and the SUSY CP problem. 
As is well-known, supersymmetry with generic CP phases is highly constrained by experimental 
limits on the neutron and electron Electric Dipole Moments (EDM). In our parameter space, 
the sfermions have multi-TeV masses and this alleviates the SUSY CP problem relative to 
realizations of SUSY with lighter sparticles. However, as we will see, constraints from EDMs are 
still relevant. 

First we consider the neutron EDM. Because the leading contribution to the neutron EDM 
is generated by the same diagram as the left of figure [2] (except with an external photon and the 
bottom quarks switched to down quarks), there is a close relationship between the size of the 
neutron EDM and the gluino contribution to the finite t/b threshold correction, 5 9 b . We note that 
the phase from the A-term is less constrained since diagrams analogous to the right diagram 
of figure [2] are suppressed by small CKM elements. For approximately degenerate squarks we 
can write the relation between the gluino contribution to the finite threshold and the neutron 
EDM EH as 




di = yArgOiMs) -^^tan/3 | -f ( -± ) . (5.5) 



i 

Since the neutron EDM is bounded by d n < 2.9 x 10~ 26 e cm [29J the bound on the CP violating 
phase goes as: 



5 



For tan j3 = 20, the LHC limit on the heavy Higgs is itla > 300 GeV [HS], which is typically weaker than the 



LHCb limit onB^ , as can be seen to the right of figure 
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where the last factor corrects for sbottom-sdown splitting. We see that a generic phase is allowed 
for multi-TeV squarks and 5( m ~ 0.1, as is required for successful b/r unification. 

The experimental limit on the electron EDM is stronger, d e < 1.05 x 10~ 27 e cm [29], and is 
generally more constraining than the neutron EDM. The leading contribution to the electron 
EDM comes from a loop mediated by charged Higgsinos and winos [67], which lead to a bound 
on Arg(/xM 2 ), 

A rg („M 2 ) < 0.02 (57) 

where, for simplicity, we have fixed /i = M 2 = 1 TeV. We see that for tan (3 ~ 50, as favored for 
t — b — t unification, the bound is stringent. The bound is alleviated for smaller values of tan/3 
(recall from above that b — r unification works for tan /3 > 10 weakening the bound on the phase 
to ~ 0.1), or with a heavier sneutrino. 



5.2 Dark Matter 

We have seen that requiring Yukawa unification along with \i < 1 TeV, in order to avoid overclo- 
sure by the thermal relic density of a neutralino LSP, places an upper bound on the superpartner 
mass scale. In this section we consider the allowed parameter space and the prospects for direct 
detection of the resulting WIMP dark matter. We shall assume that the relic abundance of the 
LSP is given by thermal freeze-out. While the parameter space of interest for Yukawa unification 
does contain single-species WIMP dark matter, we allow for the more general possibility that 
the LSP makes up only a subdominant component of the dark matter. Mult i- component dark 
matter can be obtained, for example, by environmental selection, which requires large, roughly 
equal relic densities for both WIMPs and axions due to dangerous boundaries in both the dark 
matter relic density and the vacuum misalignment angle [7J. Throughout this section we use a 
recent lattice value for the strange quark content of the nucleon [2H], f s — 0.069; some results 
are compared with a larger f s in appendix O 



Figure 17 shows the region of parameter space with D,^ < Q c d m , computed using Mi- 
crOMEGAs [68, 69J, shaded in orange, with contours of relic density normalized to Q c d m = 0.111 
(from Ref. [29]) in black. To be definite we once again fix Mi )2 using the gaugino unification re- 
lation, 6Mi fa 3M 2 ~ M3. Since the neutralino annihilation and elastic scattering cross-sections 
do not depend sensitively on M 2)3 , this assumption will not qualitatively affect our results as 
long as we have \Mi\ < |M 23 |, so that the LSP is a bino-Higgsino mixture. The relic density 
contours show that the neutralino LSP makes up an 0(1) fraction of dark matter in a large 
portion of the allowed parameter space. 

As in the previous sections, the dashed blue contours show the /i-parameter that results from 
imposing (t— )b — t unification and a 125 GeV Higgs. In both panels of the figure, the lower 
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Figure 17: The region of parameter space with $1$ < ^cdm are shown shaded in orange for 
tan /3 = 50. We have taken itla = 10 TeV in order to remove the potential effects of the heavy 
Higgs funnel and imposed gaugino unification, 6Mi = 3M 2 = M 3 . The dashed blue contours 
show the value of the /i-parameter required to achieve b — r unification, while the black contours 
show the thermal relic density of the LSP normalized to the relic abundance of dark matter, 
/Q c dm- In calculating the relic abundance, we have taken all scalar masses degenerate at mj 
and fixed A t — Af, — A T to the value required to obtain a 125 GeV Higgs. All other A-terms 
are set to 0. The darker orange shading depicts the predicted reach of the upcoming XENON1T 
direct detection experiment. 



boundary of the large orange regions corresponds to the "well-tempered" neutralino [27] . Moving 
to larger |M 3 |, the LSP becomes dominantly Higgsino-like, and so the = Q c dm boundary 
tends to follow the /i ~ 1 TeV contour, as expected for pure Higgsino dark matter. As the 
stop mass decreases, however, Higgsino dark matter is required to be heavier on account of the 
increasing importance of stop- and sbottom-mediated annihilation channels whose amplitudes 
scale as y\. If the stop is lighter still, m- t > , then it may coannihilate with the LSP in a 
small, fine-tuned region, resulting in an even heavier neutralino. 

The LSP may also be dominantly bino-like, in the case where it annihilates resonantly through 
an s— channel mediator. The Higgs and Z poles, respectively, are depicted in the left side of 
figure 17 in the upper and lower islands below the large orange region. Due to the assumption of 
gaugino unification, these regions contain a light gluino that is ruled out by LHC data; however, 
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Figure 18: Contours of direct detection cross-section are shown in purple, with the orange 



shading and blue contours as in figure 17 In terms of the elastic scattering cross-section off of 
protons, a p , we define the direct detection cross-section to be scaled by the LSP relic density, 

&DD = ®v (^jVi/^cdm) • 

relaxing that assumption would allow for a heavier gluino. 

The allowed parameter space will be further constrained by upcoming direct detection exper- 
iments, such as the planned ton-scale liquid xenon detector [70] , whose predicted reach is shaded 



in figures 17 and 18 in dark orange. However, using the lattice prediction for the strange quark 
content of the nucleon [28J, there is no limit from current direct detection experiments. Once 
again, we have used MicrOMEGAs to compute the direct detection cross-section. Contours of 
the direct detection cross-section scaled by relic abundance, a DD = a p /Sicdm) are shown 



in purple in figure 18, in which o p is the cross-section for elastic scattering of the LSP off of a 
proton. 



We see from figure [18] that XENON1T should cover the entire well-tempered region, up to 
and including = 1 TeV, at which point the LSP becomes dominantly Higgsino-like. For 
more detail on this portion of parameter space, including the expected reach of the upcoming 
XENON100 release, see appendix [C] In this region the elastic scattering is dominated by Higgs 
exchange. Increasing \M\\ reduces the Higgs-exchange cross-section, since the LSP coupling 
to the Higgs is proportional to the bino content of the neutralino. For pure Higgsino dark 
matter, the leading direct detection diagrams appear at one loop with a naive size of order 
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10~ 46 cm 2 [TjTJ [72]; however, a few groups [7JH El] find that an accidental cancellation between 
one and two- loop diagrams leads to a surprisingly small cross-section, a p < 1CT 48 cm 2 , so that 
the tree-level Higgs exchange diagrams still dominate the one and two loop contributions over 
the entire parameter range of interest. 

At sufficiently light squark masses, squark exchange diagrams with amplitudes scaling like 
yl can become important. These diagrams destructively interfere with the Higgs exchange 
contributions, and at large enough bino mass, they can dominate the elastic scattering cross- 
section. This effect accounts for the tendency of the direct detection cross-section to decrease 
and then increase again as the gaugino mass is increased with mf ~ 1 — 2 TeV. Here the 
direct detection cross-section is sufficiently small that this region of parameter space will not be 
probed by currently-planned experiments; however, UV considerations prefer a lighter gaugino 
mass. RG running to low energies with a heavy gluino pulls up on the squark masses, while 
Yukawa unification and dark matter requirements bound the squark mass to be below ~ 8.5 TeV 
or ~ 2.5 TeV, depending on the sign of X t . 

Finally, we briefly comment on the effect of relaxing the assumption of gaugino unification. If 
Mi is made larger relative to |M 3 |, then the well-tempered region moves to smaller gluino masses 
relative to the parameter space shown in figures 17 and 18, increasing the size of the parameter 
space that contains Yukawa unification without overclosure. In this case LHC searches for the 
gluino may probe some or all of the well-tempered contour. Furthermore, due to the shape of 
the /i contours, the well-tempered region would extend to lighter squark masses. Deforming in 
the other direction, in which M\ is taken lighter than M3/6, the well-tempered contour moves to 
larger gluino masses. This decreases the size of the parameter space, perhaps removing entirely 
the region with X t > 0. However, in this case direct detection experiments are able to probe 
a larger fraction of the overall space. Thus we consider the direct detection prospects in the 
interesting regions of parameter space to be promising with the next generation of experiments. 
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A Yukawa Unification with Extra Charged Matter 



Throughout this paper, we have restricted to the minimal field content of the MSSM and assumed 
that there is a desert between the superpartner mass scale and the scale of gauge coupling 
unification. One might wonder how sensitive our results are to this assumption, since the presence 
of new states will modify the beta functions of the gauge couplings and Yukawas. One possible 
addition to the field content that preserves the success of gauge coupling unification is the 
addition of complete GUT multiplets. In the context of an 577(5) GUT, for example, there may 
be some number N of extra 5 + 5's at an intermediate scale M mess . Extra charged messengers 
are typically included in models where supersymmetry is broken at a low-scale and mediated 
to the SM sector by gauge interactions. Our main focus in this paper has been non-gravitino 
DM and not low-scale SUSY breaking. Still, it is interesting to consider Yukawa unification in 
models with gauge mediation, and even if SUSY is broken at a high-scale, there may happen to 
be extra charged states present at intermediate scales. 

5 + 5 messengers 
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Figure 19: The required supersymmetric thresholds for b — r unification in the presence of extra 
charged matter. We fix tan /3 = 50 and show the values of 5^ versus ^3/^ necessary for precision 
b — r unification with TV = 0, 1, 2, 3 extra 5 + 5 with masses M mess = 100 TeV. The dashed lines 
indicate points with perfect unification, e = 0, while the shaded regions denote |e| < 0.02. The 
blue contour indicates 5b = 6[ m , assuming degenerate superpartners except for \i = 500 GeV. 
The distances between the blue contour and the dashed lines indicate the necessary values of 
the finite threshold, 5^ m , and we find that a larger finite threshold correction is required in the 
presence of extra charged states. 



In figure 19 we show the necessary value of the 5b and 63/ 5t thresholds, for precision b — r 
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unification, including N = 0,1,2,3 extra 5 + 5's at M mess = 100 TeV. We choose tan/3 = 50, 
as appropriate for t — b — r unification, and for simplicity we set 5i t 2,S T = and 5 t = 5 l ° 9 , as 
in figures [3] and [I| The blue curve indicates the value of 5 1 ^ 9 that results assuming degenerate 
superpartners except for fi = 500 GeV, as in the figures of section [2j We see that a significantly 
larger value of the finite bottom threshold is required for b — r unification in the presence of extra 
charged states. This effect is easy to understand: the new fields increase the gauge coupling beta 
functions, leading to a larger value of that pulls down harder on the bottom Yukawa, making 
the precision of b — r unification worse, before the finite bottom threshold is included. If one 
imposes bino-Higgsino DM with [i < 1 TeV, then the larger finite bottom threshold implies that 
the stop and sbottom must be even lighter than the minimal case with no extra charged states. 



B Flavor Bases 



We define here the notation we use for flavor in section 5^ We determined the flavor-dependent 
factors entering the amplitudes for b — > S7 and B s — > using the basis where the fermions 

are in mass eigenstate and the squark soft mass matrices are diagonal. In the limit of small 
left-right mixing, this corresponds to the mass basis for the squarks and, as we review below, all 
flavor violation can be parameterized by the CKM matrix and the gluino vertices [57J. In the 
following, we specialize to the (s) quark sector, which is relevant for 5-meson decays, and will 
not consider lepton flavor violation, which can be treated analogously. 

As usual, the quarks are rotated from flavor to mass basis using unitary matrices L u ^ and 

Ru,di 

u — > L u u, d —7- L d d 

u c ^R u u c , d c ^R d d c , (B.l) 

where V = L^La is the CKM matrix. Meanwhile, we rotate the squarks to the basis where the 
soft masses, ttiq, m^,m|,, are diagonal in flavor space, 

u — > L u u, d — )■ Ldd 

u c ^R u u c , d c -*R d d c , (B.2) 

where L u = L d , R u , R-d are unitary matrices. The basis with diagonal squark soft masses is 
the mass eigenstate basis in the limit of small left-right mixing. Note that left-right mixing is 
suppressed by the EW scale over the superparticle mass scale, v/rh, and therefore this basis 
is approximately equivalent to the mass eigenstate basis in the parameter space relevant for 
Yukawa unification, rh ~ 1 — 10 TeV. 

In general, different rotations are needed to bring the quarks and the squarks to mass eigen- 
state, and SUSY flavor violation arises due to this misalignment. After applying the above 
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rotations, flavor violation enters the vertices involving a quark, squark, and gaugino or Higgsino. 
Flavor violation is encoded by the following unitary matrices, 



W a = L\L a W aC = R\R a a = u,d. (B.3) 

Note that W u and W d are related by a CKM rotation: W u] W d = L\JL u V d L d = L\L d = V. The 
gluino vertices for u, u c , d, d c are simply proportional to the corresponding W a matrices. The 



charged Higgsino vertex, which also appears in section 5.1, takes the form, 



Qy u u c H u D ( dW* d V T ^-u c + dV T ^-W^uA #+ (B.4) 

where m u is the diagonal up-type quark mass matrix. The same expression describes the if7 
couplings, with the replacements u <-> d and sin/3 — > cos/3, V T —> V*. 

Much of the literature on SUSY flavor uses a basis known as the SUSY-CKM basis, where the 
entire quark superfields are rotated into the fermionic mass basis by the L U}d and R Utd matrices. 
In this basis, the squark soft masses are non-diagonal. In the limit of small flavor violation, it is 
useful to work in the mass-insertion approximation, 

T si 



M 2 d =m 2 d (t + 5 d ) 5 d = [ 6 d RL 6 d RR ) , (B.5) 



and similarly for the up-type squarks. The SUSY-CKM basis is related to the basis described 
above, with diagonal squark soft masses, by applying rotations with the W a matrices, which 
diagonalize the 5 d L and 5 d R matrices, 

diag(mj) = + W d b L d L W dV ) diag(m^) = fh 2 D {t + W dC 5% R W dcj ), (B.6) 



where the unit matrix in equation B.5 results when itlq = mp. 



C Experimental Status of The Well- Tempered Neutralino 



Here we consider in detail the well-tempered subset of the parameter space of section 5.2 



in 



which the LSP is a bino-Higgsino mixture with the correct thermal relic abundance to be all of 



dark matter. This corresponds to a portion of the boundary of the region in figures 17 and 18 
For simplicity we shall decouple all scalar superpartners, in which case the dominant elastic 
scattering diagram comes from Higgs exchange. In this case, the elastic scattering cross-section 
is proportional to |/| 2 , where 

/4 + '£/« (ci) 

9 9 / — m N 

q=u,d,s 
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MicrOMEGAs default (f s = 0.26) 



lattice (f s = 0.069) 




Figure 20: The well-tempered parameter space in the (/i,Mi) plane, with tan = 50 and 
the strange quark content of the proton taken such that / = 0.468 and 0.319 for the left and 
right sides of the figure, respectively. The green band corresponds to £1$ = 0.111 ± 0.018, 
the 3a range for the relic abundance of dark matter, while the black dashed contours show the 
elastic scattering cross-section of the neutralino against a proton. Note that, in contrast with 



the direct detection cross-sections plotted in figure 18, a p is not scaled by relic abundance. The 
shaded blue region and dashed blue contour correspond to the current limits and projected reach 
of XENON100, while the red contour depicts the expected reach of the upcoming XENON1T 
experiment. Here we decouple all scalar superpartners for simplicity, taking rh = = 10 TeV. 

Thus the ability of current and future direct detection experiments to probe the well-tempered 
region depends sensitively on the strange quark content of the nucleon, f s . Following the treat- 



ment in [75J, figure 20 shows the well-tempered parameter space in the (ju, Mi) plane, with the 
3cr range corresponding to the cosmological abundance of dark matter shaded in green. In addi- 
tion to the MicrOMEGAs default value of / s , which gives / = 0.468 [69], we show the effect of 
taking the value of f s suggested by more recent lattice calculations, which gives / = 0.319 [28J, 
considerably reducing the reach of direct detection experiments. Furthermore, in addition to 
the current XENON100 limit [76], shaded in blue, we show the predicted reach for both the im- 
minent XENON100 data release, the dashed blue contour, as well as the upcoming XENON1T 
experiment [70], drawn in red. 

Given the larger value of f s , XENON100 excludes LSP masses between m h <m t and 
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will soon be able to constrain nearly the entire well-tempered parameter space, except for very 
light neutralinos and those which annihilate resonantly through a Z or Higgs boson. However, 
if the lattice calculations are to be believed, then XENON100 only currently rules out a small 
region around ~ mw, and it will only be able to probe dark matter masses up to ~ 600 GeV. 
In either case, however, XENON1T will probe the entire well-tempered parameter space up to 
H ~ 1 TeV, at which point the LSP becomes dominantly Higgsino-like. Note that these bounds 
assume the best-case scenario: elastic scattering diagrams involving squark exchange interfere 
destructively with the Higgs exchange diagram, so that the bounds weaken as the squarks become 
lighter. 
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